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Abstract 

The possibility of determining cosmological parameters on the basis of a wide set of observational data 
including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the 
distribution of Ly-a clouds and CMB temperature fluctuations is analyzed. Using a minimization 
method, assuming JIa + il„i = 1 and no contribution from gravity waves, we found that a tilted AMDM 
model with one sort of massive neutrinos and the parameters n = 1.12 ± 0.10, il,™ — 0.41 ± 0.11 
(f^A 0.59 ± 0.11), n^drn = 0.31 ± 0.15, = 0.059 ± 0.028, VL^ = 0.039 ± 0.014 and h = 0.70 ± 0.12 
(standard errors) matches observational data best. U,^ is higher for more species of massive neutrinos, 
^ 0.1 for two and ^ 0.13 for three species. Vim raises by ~ 0.08 and ~ 0.15 respectively. The la 
(68.3%) confidence hmits on each cosmological parameter, which are obtained by marginalizing over 
the other parameters, are 0.82 < n < 1.39, 0.19 < < 1 (0 < Oa < 0.81), < fi^ < 0.17, 
0.021 < ilfc < 0.13 and 0.38 < h < 0.85. Varying only a subset of parameters and fixing the others 
shows also that the observational data set used here rules out pure CDM models with h > 0.5, scale 
invariant primordial power spectrum, zero cosmological constant and spatial curvature at a very high 
confidence level, > 99.99%. The corresponding class of MDM models are ruled out at ~ 95% C.L. It 
is notable also that this data set determines the amplitude of scalar fluctuations approximately at the 
same level as COBE four-year data. It indicates that a possible tensor component in the COBE data 
cannot be very substantial. 



1 Introduction 

The last years of the past century are marked by huge efforts of the community of astronomers, physicists 
and astrophysicists devoted to determine the most fundamental parameters of our Universe, the cosmological 
parameters. The most important among them are the mass densities of baryons (in units of the critical 
density) and of cold dark matter ilcdm, the neutrino rest masses and their total density fl^, the value of 
cosmological term A (or Ha), the Hubble constant Hq, the spatial curvature parameter Qk and the slopes n and 
amplitudes A of the primordial power spectra of scalar and tensor fluctuations. 

The primordial ratio of the number of deuterium to hydrogen nuclei (D /H) created in Big Bang nucleosyn- 
thesis is the most sensitive measure of the cosmological density of baryons fi^. Quasar absorption systems give 
definite measurements of the primordial deuterium and the most accurate value of baryon density obtained 
recently in this way is flth'^ = 0.019 ± 0.0024 0). 

The measurements of the neutrino rest mass is not so certain, unfortunately. Up-to-day we have only some 
indications for the range where it may be found. The oscillations of solar and atmospheric neutrinos registered 
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by the SuperKamiokande experiment show that the difference of rest masses between r— and ^-neutrinos is 
0.02 < Arrirfj. < O.OSeT^ 00. This also provides a lower limit for the neutrino mass, rrii, > \Am\ and 
does not exclude models with cosmologically significant values ^ 1 — 20eV. Therefore, at least two species of 
neutrinos can have approximately equal masses in this range. Some versions of elementary particle theories 
predict irii^^ « m^^ « 2.5eV and TOi^^ ~ rriiy^ ~ 10~^eV, where i^e, Vt, and Vs denote the electron, r— , /i— 
and sterile neutrinos accordingly (e.g. &). The strongest upper limit for the neutrino mass comes from the 
observed large scale structure of our Universe: ^■mi,./94eV < Q.?>h?. Since observations give for the Hubble 
parameter an upper limit of h = 0.8 one gets m,^. < 18eV. It is interesting to note that this upper limit 
coincides roughly with the upper limit for the electron neutrino mass obtained from the supernova explosion 
SN1987A and tritium /3-decay experiments. 

Important conclusions about measurements of matter density Qm (= + ^cdm + ^u) come from the Su- 
pernova Cosmology Project and the High-z Supernova Search. In particular, the relation of observed brightness 
vs. redshift for SNela shows that distant supernovae are fainter than expected for a decelerating Universe, and, 
thus, more distant. This can be interpreted as an accelerated expansion rate, or JIa > 0. The best-fit value 
is f^A = |f2m + I i 0-1 (Ic error) and flm = 1 models are ruled out at the 8a level 0. For a flat Universe 
n^ + ^A^l {^k = 0) the best-fit values are = 0.25 ± 0.1 and f^A = 0.75 ± 0.1 BBS. 

An upper limit of f^A < 0.7 (95% C.L.) follows from gravitational Icnsing statisticffl B, just consistent with 
distant supernovas results. 

Strong evidence against an open Universe can be derived from recent measurements of the position of the first 
acoustic peak in the cosmic microwave background (CMB) power spectrum by the Boomerang experiment^ . 
The la range for the curvature parameter derived from this experiment is —0.25 < fi^ < 0.15 l3' and the mean 
value is close to the flat Universe, fi^ w 0. 

Currently there are a few completely independent and broad routes to the determination of the Hubble 
constant Hq. The direct experiments can be divided into three groups: the gravitational lens time delay 
methods, the Sunyaev-Zel'dovich method for clusters and extra-galactic distance measurements. Almost all 
observations yield values of Hq in the range 50-80 km/sec/Mpc. 

Other independent methods for the determination of cosmological parameters are based on large scale 
structure (LSS) observations. Their advantage is that all parameters mentioned above can be determined 
together because the form and amplitude of the power spectrum of density fluctuations are rather sensitive 
to all of them. Their disadvantage is that they are model dependent. This approach has been carried out in 
several papers (e.g. IBB 1011) and references therein) and it is also the goal of this paper. The papers 
on this subject differ by the number of parameters and the set of observational data included into the analysis. 
In this paper a total of 23 measurements from sub-galaxy scales (Ly-a clouds) over cluster scales up to horizon 
scale (CMB quadrupole) is used to determine eight cosmological parameters, namely the tilt of the primordial 
spectrum n, the densities of cold dark matter ficdm, hot dark matter ilj^, baryons fif, and cosmological constant 
r^A, the number of massive neutrino species iV^, the Hubble parameter Hq and, in addition, the bias parameter 
bci for rich clusters of galaxies. We restrict ourselves to the analysis of spatially flat cosmological models with 
Ha + — 1 {^k = 0) and to an inflationary scenario without tensor mode. We also neglect the effect of a 
possible early reionization which could reduce the amplitude of the first acoustic peak in the CMB anisotropy 
spectrum. 

In comparison to the companion paperB the influence of the uncertainties in the normalization of the scalar 
mode amplitude caused by experimental errors on determination of cosmological parameters is also taken into 
account here. 
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2 The experimental data set and our methods 

We use the power spectrum of Abell-ACO clusters measured in the range 0.03 < k < 0.2ft,/Mpc, as 

observational input. Its amplitude and slope at lower and larger scales are quite sensitive to baryon content 
rif,, Hubble constant h, neutrino mass nii, and number of species of massive neutrinos N^, 0^ The total 
number of Abell-ACO data points with their errors used for minimization is 13, but not all of these points can 
be considered as independent measurements. Since we can accurately fit the power spectrum by an analytic 
expression depending on three parameters only (the amplitude at large scales, the slope at small scales and the 
scale of the bend); we assign to the power spectrum 3 effective degrees of freedom. 

The second observational data set which we use are the position and amplitude of the first acoustic peak 
derived from the data on the angular power spectrum of CMB temperature fluctuations. To determine the 
position and amplitude of the first acoustic peak we use a 6-th order polynomial fit to the data set on CMB 
temperature anisotropy, accumulated in Table 2 of our accompanied paper 51 data points in total. The 
amplitude Ap and position £p of first acoustic peak determined from this fit are 79.6 ± 16.5/iK and 253 ± 70 
correspondingly. The statistical errors are estimated by edges of the x^-hyper-surface in the space of polynomial 
coefRcients which corresponds to 68.3% (Icr) probability level under the assumption of Gaussian statistics. Also 
the mean weighted bandwidth of each experiment around £p is added to obtain total A£p. 

A constraint on the amplitude of the matter density fluctuation power spectrum at cluster scale can be 
derived from the cluster mass and X-ray temperature functions. It is usually formulated in terms of the 
density fluctuation in a top- hat sphere of 8h~^ Mpc radius, as, which can be easily calculated for the given 
initial power spectrum. According to the recent optical determination of the mass function of nearby galaxy 
clusters and taking into account the results from other authors (for references see B) we use the value 
asflm^^~'^''^^^'" = 0.60 ±0.08. ^From the existence of three most massive clusters of galaxies observed at z > 0.5 
a further constraint has been established by Bahcall & Fan B: S-gtlf^^ = 0.8 ± 0.1 , where a = 0.24 if JIa — 
and a = 0.29 if JIa > with Ha + = 1- 

A constraint on the amplitude of the linear power spectrum of density fluctuations in our vicinity comes 
from the study of galaxy bulk flow, the mean peculiar velocity of galaxies in sphere of radius 50/i^^Mpc around 
our position. We use the data given by KoUat & Dekel 0, V50 — (375 ± 85) km/s. 

A further essential constraint on the linear power spectrum of matter clustering at galactic and sub-galactic 
scales fc ~ (2 — 40)/i/Mpc can be obtained from the Ly-a forest of absorption lines seen in quasar spectra 0) . 
Assuming that the Ly-a forest is formed by discrete clouds of a physical extent near Jeans scale in the reionized 
inter-galactic medium at z ~ 2 — 4, Gnedin 111* has obtained a constraint on the value of the r.m.s. linear density 
fluctuations 1.6 < dpiz = 3) < 2.6 (95% C.L.) at Jeans scale for z — 3 equal to kp « 38ilm'^h/Mpc 0^. 

The procedure to recover the linear power spectrum from the Ly-a forest has been elaborated by Croft et 

al.El'. Analyzing the absorption lines in a sample of 19 QSO spectra, they have obtained the following 95% 

1/2 

C.L. constraint on the amplitude and slope of the linear power spectrum at z = 2.5 and kp — l.bilm h/Mpc 

Alikp) = klP{kp)/2Tr^ = 0.57 ± 0.26, (1) 

In addition to the power spectrum measurements we use the constraints on the value of Hubble constant 
h — 0.65 ±0.15 which is a compromise between measurements made by two groups: andB. We also employ 
the nucleosynthesis constraints on the baryon density of ilfc/i^ = 0.019 ± 0.0024 (95% C.L.) 0. 

In order to find the best fit model we must evaluate the above mentioned quantities for a given cosmological 
model. 

To this end we use the accurate analytic approximations of the MDM transfer function T(fc; z) depending 
on the parameters Q„i, ^b, ^i^, N^, h by Eisenstein & Hu 0. 
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The linear power spectrum of matter density fluctuations is 

P{k; z) = Afc"r2(fc; z)Dl{z)/Dl{Q), (3) 

where A is the normahzation constant and Di{z) is the growth factor, useful analytical approximation for which 
has been given by Carrol et alEl' . 

We normalize the spectra using the 4-year COBE data which can be expressed by the value of the density 
perturbation at the horizon crossing scale, 5h B'i-*. The normalization constant is related to 5h by 

A = 27r2<52(3000//i)3+" Mpc3+". (4) 

The Abell-ACO power spectrum is given by the matter power spectrum at z = multiplied by a linear and 
scale independent cluster biasing parameter bch which we include as a free parameter 

PA+Aco{k) = hliP{k-Q). (5) 

For a given set of parameters n, Sim, S^b, h, O^, and bd the theoretical value oi PA+Aco{kj) can now be 
calculated for each observed scale kj. Let's denote these values by yj {j = 1, 13). 

The dependence of position and amplitude of the first acoustic peak in the CMB power spectrum on cos- 
mological parameters has been investigated using the public code CMBfast by Seljak & Zaldarriagall' . As 
expected, these characteristics are independent on the hot dark matter content. We determine the values ip 
and Ap for given parameters (n, h, Qb and SIa) on a 4-dimensional grid for parameter values in between the 
grid points we determine £p and Ap by linear interpolation. We denote £p and Ap by yu and yi5 respectively. 

The theoretical values of the other experimental constraints are obtained as follows: The density fluctuation 
as is calculated according to 

1 

al^— k^P{k; 0)W^{8Mpc k/h)dk, (6) 
27r^ Jo 

with P{k;z) from Eq. (H). We set yia = crg^^™'*'^"""^"'" and yn = crgll", where a = 0.24 for f^A = and 
a — 0.29 for Qa > 0, respectively. 

The r.m.s. peculiar velocity of galaxies in a sphere of radius R — 50/i^^Mpc is 

V^Q = — k^P^''\k)e-'''^'fW^{50Mpck/h)dk, (7) 
27r^ Jo 

where P^^\k) is the density-weighted power spectrum for the velocity field 0', M^(50Mpc k/h) is a top-hat 
window function, and Rf = 12/i~^Mpc is the radius of a Gaussian filter used for smoothing of the raw data. 
For the scales considered P'-'"^{k) « {Q^-^ Hq)^ P{k; 0) / k"^ . We denote the r.m.s. peculiar velocity by yig. 

The value of the r.m.s. linear density perturbation from the formation of Ly-a clouds at redshift z and scale 
kp is given by 

al{z) = —j eP{k-z)e^-''"'^^"dk. (8) 

We set cr^(z = 3)?/i9. 

The value of Ap(fcp,z) and the slope n{z) are obtained from the linear power spectrum P{k;z) by Eq. (|l|) 
and Eq. (||) at z = 2.5 and kp — 0.008i?(z)/(l -I- z)(km/s)^^, and are denoted by y2o and y2i accordingly. 

For all tests except Gnedin's Ly-a test we use the density weighted transfer function Tcbu{k, z) from El*. For 
(jp the function Tcb{k, z) is used according to the prescription given by Gnedin 0. Note, however, that even in 
the model with maximal ^li, (~ 0.2) the difference between Tcb{k, z) and Tcbv{k, z) is less than 12% for k < kp. 

Finally, the values of flbh^ and h are denoted by ?/22 and y23 respectively. 

Under the assumption that the errors on the data points are Gaussian, the deviations of the theoretical 
values from their observational counterparts can be characterized by x^' 

23 / ~ X 2 

^2 _ / Vj - yj \ (ct^ 
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where yj and Aj/j are the experimental data and their dispersions, respectively. The set of parameters n, 
flm, ^b, h, f2y, Nu and hd are then determined by minimizing using the Levenberg-Marquardt method 
. The derivatives of the predicted values w.r.t the search parameters required by this method are calculated 
numerically using a relative step size of 10^^. 

This method has been tested and has proven to be reliable, independent on the initial values of parameters 
and it has good convergence. 



3 Results 

The determination of the parameters n, ri,„, fif,, h, Q,^, N^^ and b^i by the Levenberg-Marquardt minimization 
method is realized in the following way: we vary the set of parameters n, Qm, ^b, h, Q^, and bd and find the 
minimum of x^, using all observational data described in previous section. Since the N^, is discrete, we repeat 
the procedure three times for N^=l, 2, and 3. The lowest of the three minima is the minimum of for the 
complete set of free parameters. 

Table 1: Cosmological parameters determined for the tilted AMDM model with one, two and three species of 
massive neutrinos. 



Arnin 


n 








h 


bci 


1 4.64 

2 4.82 

3 5.09 


1.12±0.09 
1.13±0.10 
1.13±0.10 


0.41±0.11 
0.49±0.13 
0.56±0.14 


0.059±0.028 
0.103±0.042 
0.132±0.053 


0.039±0.014 
0.039±0.014 
0.040±0.015 


0.70±0.12 
0.70±0.13 
0.69±0.13 


2.23±0.33 
2.33±0.36 
2.45±0.37 



The results are presented in the Table 1. The errors in the determined parameters are the square roots 
of diagonal elements of the covariance matrix of standard errors. More information about the accuracy of 
the determination of parameters and their sensitivity to the data used can be obtained from the contours of 
confidence levels presented in Fig. |l] for the tilted AMDM model with parameters from Table 1 (case = !)• 
These contours show the confidence regions which contain 68.3% (solid line), 95.4% (dashed line) and 99.73% 
(dotted line) of the total probability distribution in the two dimensional sections of the six-dimensional parameter 
space, if the probability distribution is Gaussian. Since the number of degrees of freedom is 7 they correspond 
to Ax^ =8.2, 14.3 and 21.8 respectively. The parameters not shown in a given diagram are set to their best-fit 
value. 

As one can see in Fig. la the iso-x^ surface is rather prolate from the low-flm - high-n corner to high-Om 
- low-n. This indicates a degeneracy in n — Qm parameter plane. Within the la the 'maximum likelihood 
ridge' in this plane can be approximated by the equation n^/TT^ — 0.73. A similar degeneracy is observed in 
the ill, — plane in the range < fti, < 0.17, 0.25 < flm < 0.6 (Fig.lc). The equation for the 'maximum 
likelihood ridge' or 'degeneracy equation' has here the form: il^ — 0.023 — 0.4451„i -I- 1.3f2^j . 

The important question is: which is the confidence limit of each parameter marginalized over the other ones. 
The straightforward answer is the integral of the likelihood function over the allowed range of all the other 
parameters. But for a 6-dimensional parameter space this is computationally time consuming. Therefore, we 
have estimated the la confidence limits for all parameters in the following way. By variation of all parameter 
we determine the 6-dimensional x^ surface which contains 68.3% of the total probability distribution. We 
then project the surface onto each axis of parameter space. Its shadow on the parameter axes gives us the la 
confidence limits on cosmological parameters. For the best AMDM model with one sort of massive neutrinos 
the la confidence limits on parameters obtained in this way are presented in Table 0. 
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Figure 1: Likelihood contours (solid line - 68.3%, dashed - 95.4%, dotted - 99.73%) of the tilted AMDM model 
with Nv = 1 and parameters from Table 1 {N^ = 1) in the different planes of n — Sim — Cli, — flb — h space. The 
parameters not shown in a given diagram are set to their best fit value. 
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Table 2: The best fit values of all the parameters with errors obtain by maximizing the (Gaussian) 68% 
confidence contours over all other parameters. 



parameter 


central value and errors 


n 
bci 


41+0-59 
U.UO_o.06 

U.UO»_Q 

^ -p,+0.15(+0.31) 
'^-0.32 

1 19+0.27 
-■^•^^-0.30 

00+1.3 



*^ - the upper limit iS|-abtained by including the lower limit on the age of the Universe due to the age of oldest 
stars, to > 13.2 ± 3.0 IliP. The value obtained without this constraint is given in parenthesis. 



It must be noted that, using the observational data described above, the upper la edge for h is equal 1.08 
when we marginalized over all other parameters. But this contradicts the age of the oldest globular clusters 
to — 13.2 ± 3.00. Thus we have included this value into the marginalization procedure for the upper limit of 
h. We then have 8 degrees of freedom (24 data points) and the 6-dimensional surface which contains 68.3% 
of the probability is confined by the value 13.95. We did not use the age of oldest globular cluster for searching 
of best fit parameters in general case because it is only a lower limit for age of the Universe. 

The errors given in Table ^ represent 68% likelihood, of course, only when the probability distribution is 
Gaussian. As one can see from Fig.l (all panels without degeneracy) the ellipticity of the likelihood contours 
in most of planes is close to what is expected from a Gaussian distribution. This indicates that around their 
maxima the likelihood functions are close to Gaussian. However, the asymmetry of the error bars obtained, 
shows that away from the maxima this is no longer the case. Therefore, our estimates of the confidence limits 
have to be taken with a grain of salt. The errors define the range of each parameter within which the best- 
fit values obtained for the remaining parameters lead to Xmin ^ 12.84. Of course, the best-fit values of the 
remaining parameters lay within the range given in Table ^j. However clearly not every set of parameters from 
these ranges satisfies the condition, Xmm — 12.84. For example, standard CDM model {il„i — 1, h — 0.5, 
rib — 0.05, n = 1 and best-fit value of cluster biasing parameter bd — 2.17 (tis — 1-2)) has Xmin — 142 (!), 
which excludes it at very high confidence level, > 99.999%. When we use the baryon density inferred from 
nucleosynthesis {h'^flb = 0.019 {bd = 2.25, as = 1.14)) the situation does not improve much, Xmin = 112- 
Furthermore, even if we leave h as free parameter we still find Xmin — 1^ l'^) with the best-fit values 
h ~ 0.37 and bd — 3.28 {as — 0.74); this variant of CDM is ruled out by direct measurements of the Hubble 
constant. 

The standard MDM model (0„ ^ 1, h = 0.5, fib = 0.5, n = 1, fi^ = 0.2, = 1 with a best value of the 
cluster biasing parameter bd = 2.74 (erg = 0.83)) does significantly better: it has Xmin = 23.1 (99% C.L.) which 
is out of the 2a confidence contour but inside 3a. With the nucleosynthesis constraint the situation does not 
change: Xmin — 22; also if we leave h as free parameter: Xmin — 21) h — 0.48. But if, in addition, we let vary 
ri^, we obtain Xmin — 13 with best-fit values of il^, — 0.09, h = 0.43, bd — 3.2 {as = 0.73). This means that 
the model is ruled out by the data set considered in this work at ~ 70% confidence level only. But also here 
the best-fit value for h is very low. If we fix it at lower observational limit h = 0.5 then Xmm — 1^-9 (the best 
fit values are: fi^ = 0.15, bd = 2.8 (erg = 0.83)), which corresponds to a confidence level of 95% . 
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Therefore, we conclude that the observational data set used here rules out CDM models with h > 0.5, a 
scale invariant primordial power spectrum (n = 1) and Qk = f^A = at very high confidence level, > 99.99%. 
MDM models with h > 0.5, n = l and ^ = are ruled out at 95% C.L. 

One can see the model with one sort of massive neutrinos provides the best fit to the data, Xmin ~ 4.6. Note, 
however, that there arc only marginal differences in Xmin ^'^^ = 1:2,3. Therefore, with the given accuracy 
of the data we cannot conclude whether - if massive neutrinos are present at all - their number of species is 
one, two, or three. 

The number of degrees of freedom is Np = N^^p — Nps,^ — 7. The Xmin all cases is within the expected 
range, Np — \/2Nf < x^in — + V^Np for the given number of degrees of freedom. This means that the 
cosmological paradigm which has been assumed is consistent with the data. 

One important question is how each point of the data influences our result. To estimate this we have 
excluded some data points from the searching procedure. Excluding any part of observable data results only 
in a change of the best-fit values of n, Qm and h within the range of their corresponding standard errors. This 
indicates that the data are mutually in agreement, implying consistent cosmological parameters (within the still 
considerable error bars). The small scale constraints, the Ly-a tests reduce the hot dark matter content from 
ri,y ~ 0.22 to ~ 0.075. The cs-tests further reduce fijy to ~ 0.06. Including of the Abell-ACO power spectrum 
in the search procedure, tends to enhance f2^ slightly. The most crucial test for the baryon content is of course 
the nucleosynthesis constraint. Its ~ 6% — Icr-accuracy safely keeps h^ilb near its median value 0.019. The 
parameter fib in turn is only known to 36% accuracy due to the large errors of other experimental data used 
here, especially of the Hubble constant. The accuracy of h (~ 17%) is better than the one assumed from direct 
measurements, ~ 23%. Summarizing, we conclude that all data points used here are important for searching 
the best-fit cosmological parameters and do not contradict each other. 

Up to this point we ignored the uncertainties in the COBE normalization. Indeed, the statistical uncertainty 
of the fit to the four- year COBE data, 6^, is 7% (Icr) B and we want to study how this uncertainty influences 
the accuracy of cosmological parameters which we determine? 

Varying Sh in the la range we found that the best-fit values of all parameters except Q^, do not vary by 
more than 2% from the values presented in Table 1. Only U^, varies in a range of 12% . These uncertainties are 
significantly smaller than the standard errors given in Table 1 and neglecting them is thus justified. The normal- 
ization constant A, which for best model with one species of massive neutrinos, A = 4.68 • 10^(/i^^ Mpc )'^^", 
varies in a range of 8%, and not 14% as one might expect from (^) at the first sight. The reason is that vari- 
ation of 6h is somewhat compensated by correlated variation of n and h. Moreover, if we disregard the COBE 
normalization and treat the normalization constant A as a free parameter to be determined like the others, its 
best-fit value becomes 4.82- 10'^{h~^ Mpc )^+" (for — 1), consistent with COBE normalization. The best-fit 
values of the other parameters correspondingly do not vary substantially: n = 1.09, fi™ = 0.40, fli, ~ 0.052, 
rib — 0.041, h — 0.68 and bd — 2.19 (this is less than 5% except for il^, which is reduced by 12%). This 
implies that determinations of the amplitude of scalar fiuctuations by the COBE measurement of the large 
scale CMB anisotropics and by large scale structure data at much smaller scales are in good agreement. It also 
indicates that a possible tensor component in the COBE data cannot be very substantial. 

4 Conclusions 

We summarize, that the observational data of the LSS of the Universe considered here can be explained by 
a tilted AMDM infiationary model without tensor mode. The best fit parameters are: n — 1.12 ± 0.09, 
= 0.41±0.11, = 0.06±0.028, fib = 0.039±0.014and h = 0.70±0.12. AU predictions of measurements are 
close to the experimental values given above and within the error bars of the data. The CDM density parameter 
is ^cdm = 0.31 ± 0.12 and is moderate, = 0.59 ± 0.11. The neutrino matter density corresponds to a 
neutrino mass m,y — 94il^/i^ « 2.7 ± 1.2 eV. The value of the Hubble constant is close to the measurements 
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by Madore et al.O. The age of the Universe for this model equals 12.3 Gyrs which is in good agreement with 
the age of the oldest objects in our galaxy 0. The spectral index coincides with the COBE prediction. The 
relation between the matter density Vim and the cosmological constant J^a agrees well with the independent 
measurements of cosmic deceleration and global curvature based on the SNIa observation. 

The 1(7 (68.3%) confidence limits on each cosmological parameter, obtained by marginalizing over the other 
parameters, are 0.82 < n < 1.39, 0.19 < f7„ < 1, < f^A < 0.81, < < 0.17, 0.021 < fib < 0.13 and 
0.38 <h< 0.85. 

The observational data set used here rules out the CDM models with h > 0.5, scale invariant primordial 
power spectrum n = 1 and Oa = fife = at very high confidence level, > 99.99%. Also pure MDM models are 
ruled out at - 95% C.L. 

It is remarkable also that this data set determines the value of normalization constant for scalar fluctuations 
which approximately equals the value deduced from COBE four-year data. It indicates that a possible tensor 
component in the COBE data cannot be very substantial. 

The coincidence of the values of cosmological parameters obtained by different methods indicates that a wide 
set of cosmological measurements are correct and that their theoretical interpretation is consistent. However, 
we must also note that the accuracy of present observational data on the large scale structure of the Universe 
is still insufficient to determine a set of cosmological parameters with high accuracy. 
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